^ 






(N 



Quantum Energy Teleportation with a Linear Harmonic Chain 



Yasusada Nambio 

Department of Physics, Graduate School of Science, 

Nagoya University, Chikusa, Nagoya 464-8602, Japan 



I 



Masahiro Hott 

Q ! Department of Physics, Faculty of Science, 

Tohoku University, Sendai 980-8578, Japan 
O 
Q ■ (Dated: October 14, 2010) 



Abstract 

A protocol of quantum energy teleportation is proposed for a one-dimensional harmonic chain. 
A coherent-state positive operator-valued measure (POVM) measurement is performed on coupled 
oscillators of the chain in the ground state accompanied by energy infusion to the system. This 



qh' measurement consumes a part of the ground state entanglement. Depending on the measurement 



result, a displacement operation is performed on a distant oscillator accompanied by energy ex- 



>r , traction from the zero-point fluctuation of the oscillator. We find that the amount of consumed 

'NT . 

L -J I entanglement is bounded from below by a positive value that is proportional to the amount of 

• ■ teleported energy. 
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I. INTRODUCTION 



Recently, it has been reported that energy can be transported by local operations and clas- 



sical communication (LOCC) w 

causa 

tems 



ity ( for spin systems 

3B). 



hile retaining local energy conservation and without breaking 
pj-yl, for trapped ion systems |^, and for quantum field sys- 
7|). Such protocols are called quantum energy teleportation (QET) and are based 
on ground-state entanglement of many-body quantum systems. By performing a local mea- 
surement on a subsystem A of a many-body system in the ground state, information about 
the quantum fluctuation of A can be extracted. During this measurement, some amount 
of energy is infused into A as QET energy input, and the ground state entanglement gets 
partially broken. The measurement result is announced to another subsystem B of the many- 
body system by the classical protocol. Using this information, energy can be extracted from 
B by performing a local operation on B dependent on the announced measurement data. 
The root of the protocols is a correlation between the measurement information of A and 
the quantum fluctuation of B via the ground state entanglement. The information about 
A enables us to partially know what kind of zero-point fluctuation of B is realized. Thus, 
using this information, we can select a good operation on B for the energy extraction. 

In general, we are able to make a better strategy for a task by obtaining more infor- 
mation. Hence, for the QET case, it sounds plausible to imagine that more information 
about the quantum fluctuation of B is obtained by measurements of A, more energy can 
be teleported from A to B. If we consume a large amount of ground-state entanglement 
between fluctuation of A and fluctuation of B during the measurement of A, it is naturally 
expected that much information about the post-measurement state of B is included in A's 
measurement result. Therefore, there should exist some qualitative relationship between 
the breaking of entanglement by measurement and the amount of teleported energy of the 
optimal operation of B. 

While the basic protocol of QET has been investigated for several model systems, the 
analysis of QET for many-body systems has not been development yet. This is because these 
quantities are difficult to obtain analytically for many-body systems. So far, the relation 
between teleported energy and entanglement breaking has been investigated only for the 
minimal model of QET ^] which consists of two qubits. In this paper, we aim to establish 
this relation of QET using a one-dimensional harmonic chain. There are many works in- 



vestigating the feature of entanglement in harmonic chain models |8Hl2|. The conventional 
method to analyze the entanglement in a harmonic chain is to introduce two spatial regions 
containing sites of the harmonic chain and consider the bipartite entanglement between these 
regions. Previous analysis shows that the zero temperature ground state (vacuum state) of 
the harmonic chain is entangled. Thus this entanglement can become a resource for QET. 
As we are not interested in the dynamical aspect of QET in the present analysis, we assumed 
that the speed of classical communication is infinity and our protocol can be treated as the 
non-relativistic one. 

For a harmonic chain, we first consider the setting that groups A and the B consist of 
a single site. We investigate the relation between the amount of teleported energy and the 
quantum mutual information by changing the distance between A and B. Then we consider 
the setting with a block of coupled harmonic oscillators as A and its complementary set 
as B. Increasing the number of measured oscillators in A, we calculate both the entangle- 
ment breaking between A and B and the amount of energy teleported from A to B via a 



QET protocol. As a measure of entanglement, we use the logarithmic negativity ISj. An 
explicit inequality is given such that the breaking of entanglement is lower bounded by a 
positive value that is proportional to the amount of teleported energy. This shows how much 
entanglement is required to teleport energy in the harmonic chain. 

The paper is organized as follows. In Sec. II, we introduce the harmonic chain model 
and we prepare a formula for the positive oeprator-valued measure (POVM) measurement 
for the harmonic chain in Sec. III. We present our numerical result in Sec. IV and Sec. V is 
devoted to a summary. 

II. HARMONIC CHAIN MODEL 

We consider a one-dimensional harmonic chain with N sites. We assume that A^ is an 
even number. The Hamiltonian is 

1 ^ 
^=25Z(P? + ^|-"^^-^^-0' (j = l,---,Ar) (1) 

where we assume a periodic boundary condition for qj,Pj and a positive parameter a 7^ 1 
is introduced to regularize the infrared divergence which appears in the correlation matrix. 
a = 1 corresponds to the critical case (massless limit) of the harmonic chain. This criticality 



allows us to teleport energy to a distant point, decreasing the amount of energy obeying 
a power law decay with respect to the distance. As we will show, QET is possible for the 
non-critical case. However, the amount of teleported energy decreases exponentially with 
increasing distance and QET is not so effective. Thus, it is important to investigate QET 
with near-critical harmonic chains. The quantized canonical variables are 

1 ^-^ 



N-l 

-^ E (^^'^ + /^*«!v-.) e^''=^ (2) 



^^^ fc=0 



where fk,gk are mode functions 



1 -iLOkt n. - . \^^-i^ht ,,,2 



j^ = -=e-^-^\ g^ = ^^e-'--\ ujl = \-aco^Qk (3) 



C^fc? ^y 



and the creation and annihilation operators satisfy 

^kk'- 

The two point correlation functions with respect the ground state satisfying ak\g) = are 

1 ^"' 1 

{qiqj) = j^^^^o^ [(^ - J)^k] = g\t-j\ = Gij, (4) 

k=0 ^ 

^ N-l 

fc=0 
. N-l 

A;=0 
{(li) = {Pi) = 0, 

where we have introduced the matrices G and H satisfying the relation 

{GH),j = ls.r (7) 

As the ground state is a Gaussian state, it can be completely determined using these two 
point correlation functions. 

III. POVM MEASUREMENT AND ENERGY TELEPORTATION 

We briefly review the protocol of QET (Fig. [1]). In these figures, quantum fluctuations 
in the harmonic chain are schematically shown as broken lines. The dashed horizontal lines 



represent the magnitude of quantum fluctuation of the state. The protocol consists of the 
following steps: 




Ea=0 
(2) measurement of A 



Ea>0 



Eb=0 



N. Announce 

•^ the result to Bob 




Local Unitary 
operation 

Ub(X,P) 




Eb <0 



FIG. 1: The protocol of QET with a harmonic chain. (1) We prepare the ground state of a 
harmonic chain and introduce two parties A and B. The local energy of A and B is zero. (2) 
The measurement of A increases the energy of the group A. The result of the measurement is 
announced to Bob via classical communication. (3) Depending on the announced measurement 
result, Bob performs a local unitary operation on the group B, which reduces the energy of the 
group B so that it becomes negative. The dashed horizontal lines represent the magnitude of 
quantum fluctuation of the state. 

6 



(1) We prepare the ground state of the harmonic chain and define groups A and B of coupled 
oscillators in the system. A and B consist of some number of sites. The groups A and B are 
entangled due to the entanglement of the ground state quantum fiuctuation of the harmonic 
chain. This entanglement enables us to extract a positive energy via a protocol of QET. 

(2) Alice makes a local measurement of quantum fiuctuations of sites in the group A. 
Because the post-measurement state is not the ground state, it is an excited state with 
positive energy. As some amount of energy is infused through this measurement procedure, 
the energy of the group A after the measurement becomes positive provided that the ground 
state energy before the measurement is chosen to be zero. By this measurement, Alice 
obtains measurement results for sites in the group A. She announces these values to Bob 
via classical communication. 

(3) Depending on the announced measurement results. Bob performs a local unitary 
operation on sites in the group B. By choosing the operation suitably to suppress the zero- 
point fiuctuation of B, it can be shown that the local energy of the group B takes a negative 
value lower than the ground-state one. Thus we can extract positive energy from the system. 

In our numerical investigation of QET for the harmonic chain, we apply a POVM mea- 
surement [ij] to sites in the group A (see Fig. [2] and [6]). The measurement reduces the 
entanglement between the two groups. What we want to know is the relation between the 
amount of entanglement breaking due to the measurement and the extracted energy via the 
protocol of QET. 

A. The state before and after the measurement 

We assume that the state of the system before the measurement is the ground state. In 
the coordinate representation, the state is given by 

po = {q\9){9\q), q = {qir-- ,(lNf (8) 

where 



{q\g) oc exp 



exp 



N 



^ qiHijqj 



-- J2 q^i.G %qj 

We apply a coherent-state POVM measurement to 2£ + 1 sites in the group A. To define 
the POVM adopted in our analysis, let us introduce the following annihilation operator for 



the canonical variable {qn,Pn) of the site n in the group A: 



Qn 



Pr. 



1 < n<2£+l. 



2 - V2W 
We define the ground state |0„) with respect to the operator 6„ as 

and the coherent state defined by b}^ as 

|X„,P„) = e-l'^-l'/^exp (cJi) |0„), c„ = ./|x„ + ^P„, (9) 

V / V ^ v2(^ 

where X„ and P„ are real numbers and are classical amplitude of the coherent state. The 
expectation values of site variables with respect to this coherent state are 



{Xn-, Pn\<in\Xn-, Pn) = -^n + ^5 {^n-, Pn\Pn\Xn-, Pn) = -P„ + — • 

In the coordinate representation, the coherent state can be written 



(10) 



l^ra; -* n/ 



CJ\l/4 



TT 



exp 



UJ 



— Xn) +iP„(g„ — X„ 



The measurement operator for the POVM is given by 

2^+1 



where 



2^+1 N 

n|X„,P„)(X„,P„|, l<i< 2 

^ ' "' "' \ ri-, n\-, _ _ 2 

' n=l 



(11) 



(12) 



X — (Xi, ■ ■ ■ , X2^+i) , P — (Pi 
The POVM is defined by 



, ■ ■ ■ 5-^2^+1; 



2£+l 



n,(X,P) = Ml{X,P)M,{X,P) = — -i^ J] |X„,P„)(X„,P„| 

^ ^ n=l 

and satisfies the following completeness relation for the measurement operator: 



(13) 



/"( fjrfx„rfpjn,(x,p) = j. 



If we obtain the value (X, P) as the result of the POVM measurement, the state after the 
measurement can be written as 

2^+1 2^+1 

P,{X,P) = l[{X^,Pn\9){9\Xn,Pn)^l[ |X„,, P„,) (X„,, P„, | . (14) 



n=l 



n'=l 



After the measurement, the states for the measured sites become the coherent states ([9]) 
and states for other sites are computed by acting the measurement projection operator on 
the ground state. To obtain the exphcit form of the state after the POVM measurement, 
we introduce the following vector notation for the site variables: 

(Qn) = r' I , {Pn] ' 

The (2£ + l)-dimensional vectors qf^,pf^ represent variables for the sites whose values are 
measured. Then, with respect to these bases, the correlation matrices (jl]) and (jSD for site 
variables can be written 



H= \ \ , G 





where Li, Ci denote the {2i+ 1) x {2i+ 1) correlation matrices for the measured sites. Using 
this notation, the ground state projected to the coherent state (jH]) is given by 



{X,P\g) ocexp 



-qj(H,-Kj(^Le + ^) ' K^ q, + (^ KJ (^L, + ^^ \iP - uX) 



The two-point correlation functions for fluctuations Ag^ = Qi — (qi), ^Pi = Pi — (Pi) after the 
measurement are obtained by taking the expectation value with respect to the state (fT^ . 
The non-zero correlation functions are 

i,j e (measured site in A) : (AqiAqj) = -^, {ApiApj) = -Sij, 

{ApiAqi + AqAPi) = 
i,j i (measured site in A) : {Aq^Aq^) = - {M''^)i_^2e+i),j-(2e+i) ' (1^) 

(ApiApj) = (M)._(2^^i)_^._(2f^^) 



where the matrix M is defined by 



M = He-K[{Le+'^) K,. (17) 



The other two-point functions are zero. These covariances are independent of the result of 
the measurement {X, P). The probabihty distribution to obtain the measured value (X, P) 
is given by 



(^|n,(X,P)|^)= I dq,\{g\X,P)\' 

1 



oc exp 



B. The optimized energy 



We define the energy of the lattice site n as 



-^P^iLe + ^yp-lx^(ce + ^] X 



-1 



(18) 



Hr, 



a 



Pn + ln- -^<in[qn-l + Qn+l^ 



— e 



(19) 



where the constant e is chosen to satisfy 



{g\H^\g) = 0. 



(20) 



After the measurement, we obtain the values (X, P) of the measured sites. Depending on 
these measured values, we apply the following local displacement operation to sites in the 
group B which is composed of a single site in our numerical setting {ub = N/2 + £ + 1) 



f/B(X,P) = exp 



P)qn, - (0 ■ X)p.. 



riB 



(21) 



where 9 and are constant parameters to be determined later by optimization. Then, the 
transformed state is 

,2^+1 



P', 



fl[{dXr,dP^)UB{X,P)M,{X,P)\g){g\Ml{X,P)Ul{X,P). (22) 

^ n=l 



The energy of the site B with respect to this state is 

2^+1 2€+l 



Mr-1, 



1^) 



{Hb) = Y^ OjiglPjPnslg) + Yl ^Aalij (^ns - -^((ins+i + q^ 

+ 2 5Z ^J' [(9\PjPk\g) + l^ifc) dk + 'Y^f^A ialQMg) + — J 

7,fc=l i,k=l ^ ^ 



4>k 



\d^Tj+ J, ■ e + ^^tJ+ jg ■ 



(23) 
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where we have introduced 



(^p)jfc - ^\J-k\ + 2 ^J^' C^i^jk - 9\j-k\ + ^ Sjk, 



(24) 



Jr 



P J 



hi 



Jr, 



'|j-"sl' \yq j 3 ~ 9\j-nB\ 

By choosing the following values of the parameters 



f|j„(„^_i)l +5f|j_(„^+i)|) . 



U ip Jp, 






(25) 



we obtain the following minimum (optimized) value of the site 5's energy 

{HB)opt - 



2 P P P 2 '^ 1 '^1^^- 



(26) 



The negative sign of {Hb) means that a positive amount of energy is transferred from the 
chain to outside, obeying local energy conservation. 

C. Bipartite entanglement and mutual information 

To investigate the relation between the amount of extractable energy via QET and the 
entanglement breaking due to the measurement, we consider the bipartite entanglement 
between the group A and the group B. For this purpose, we use the logarithmic negativity 
En as our measure of entanglement. The logarithmic negativity provides an upper bound 
on the efficiency of distillation of entanglement 13|. 

The canonical variables of our bipartite system are 



?-{ 



qA, PA, QB, pb 



where O is a 2N x 2N matrix 



n 



( 



\ 



y [0'^fc] = «fiifc, j, A; = I,--- ,2A^ 



\ 



7 




As we are assuming a Gaussian state, the state is completely specified by the following 
covariance matrix 

V^.fc = ^fea + aO)-(0)(a)- (27) 
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The symplectic eigenvalue z/^ of the covariance matrix V satisfies the inequahty 

^n>l (28) 

which is the generahzation of the uncertainty relation. The logarithmic negativity is defined 
by the partially transposed covariance matrix V obtained by reversing the sign of party 5's 
momentum 

N 

En = -J2 ^"^ [0' log2 (2^>n)] (29) 

n=l 

where z/„ is the symplectic eigenvalue of the partially transposed covariance matrix V. If 
this quantity is positive, the bipartite system is entangled, and we can use the logarithmic 
negativity as the measure of the bipartite entanglement between groups A and B. The 
logarithmic negativity provides the sufficient condition for entanglement between A and B. 
We also consider the quantum mutual information of A and 5 as a measure of correlation 
between A and B. For the Gaussian system with the covariance matrix V, in terms of 



symplectic eigenvalues ly^, the von Neuman entropy of the system is given by [15| 

For the bipartite system with the group A and B, the mutual information is defined by 

SM = S{A) + SiB)-SiA + B). (31) 

This quantity represents the total correlations including the quantum and the classical cor- 
relation between the group A and B [16]. A large value of Sm implies that A and B are 
strongly correlated and we are able to obtain much information about B just by measuring 
A. Thus, this is expected to assure a large amount of teleported energy. 

We numerically calculate the logarithmic negativity and the mutual information. The 
covariance matrix of the total system before the measurement is 



Vn 



Uq • • • (^0 



jj(N,l) rj{N,N) 

\Uq • • • (^0 / 



{ApiApj) \ /i|,_,| 



The covariance matrix after the measurement is 



T/M n 

F^ = I ^ I , (33) 

Vab. 
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where V^ represents the covariance matrix for the measured sites and is given by the 
2(2£ + 1) X 2(2£ + 1) diagonal matrix 



/ 



V, 



M 



Ui 



\ 



\ 



Ui 



f/l 




(34) 



and Vab represents the covariance matrix for the un-measured sites, 

/ jj{l,l) ___ jjil,N-2d-l) \ 



Vab 



Tj{N~2d-l,l) Tj{N-2d-l,N-2d~l) 



'1 



, [/(*'■'■) 



(M-1) 







4 ^""" '^<o 

(M), 



(35) 



«j/ 



By use of these expressions for the covariance matrices, the logarithmic negativity and the 
mutual information before and after the measurement can be obtained numerically. 

IV. NUMERICAL RESULT 

We numerically calculated the decrease of the logarithmic negativity due to the measure- 
ment 

Ai^jv = -Eat (before measurement) — -Ejv (after measurement) 

and the optimized energy of the site B obtained via the protocol of QET. We consider 
the two different settings of numerical calculations to establish the relation between the 
entanglement breaking and the amount of teleported energy. As the parameter a of the 
harmonic chain, we used the following four values: 



ai 



0.90, ^2 = 0.95, ^3 = 0.99, ^4 = 1 - 10"^ (critical). 



(36) 



We have also investigated harmonic chains with smaller values of a, but the amount of 
teleported energy is reduced and QET is not effective for small values of a. We adopt these 
four values of a as typical cases of our numerical calculation. The values ai, 0:2, cts represent 
typical examples of the non-critical case and 0:4 represents the critical case with a small 
cutoff parameter 10^^ which is necessary to evaluate correlation matrices numerically. As 
we will see later, the numerically obtained quantities converge to the value with 0:4 as the 
value of a approaches 1. 
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A. Setting 1 



A^ 




FIG. 2: The setting 1 of our numerical calculation. We change the separation d between sites A 
and B. 



As the first setting of our numerical calculation, group A and group B consist of a single 
site (Fig. 12]). We fix the total site number A^ = 100. We apply the POVM measurement to 
the site A. By changing the separation d between A and B, we observe how the entanglement 
between the two groups and the optimized energy (12 6 p of the site B change. 

0.5 r 



< 




0.94 0.96 
a 



1.00 



FIG. 3: Left panel: d dependence of the decrease of the logarithmic negativity due to the measure- 
ment for a4 = 1 — 10~^. The logarithmic negativity before and after the measurement is zero for 
d > 1 and this behavior is the same for the other values of a. Right panel: a dependence of the 
AEn for d = 0. 



Figure |3] shows AE^ as a function of the separation between A and B. The bipartite 
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system composed of A and B is entangled only for d = before the measurement. After 
the measurement, A and B become separable even for d = 0. This behavior is the same 
for other values of a. As a approaches 1, the entanglement consumed via measurement 
increases. Our numerical calculation indicates that we have no entanglement breaking for 
d = 1, 2, 3, ■ ■ • . However, AEj^j = for d ^ does not mean sites A and B do not have the 
correlation necessary to establish the QET protocol. 




FIG. 4: d dependence of the optimized energy of the site B. For the critical case a = 04, this 
quantity behaves as \{Hb)\ ~ 2 x 10^^d~^'^ for large d. 

Indeed, as Fig. H] shows, we have non-zero (negative) optimized energy of the site B even 
for AEi^ = {d > 1). The amount of optimized energy of site B decreases as d increases. 
For the critical case a^, the optimized energy for large d behaves as 



(H 



B 



2 X 10"^rf" 



-3.6 



(37) 



To understand the behavior of the optimized energy from the view point of entanglement 
breaking, we investigated the d dependence of the mutual information f l?!]) (Fig. [S]). 
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FIG. 5: d dependence of the decrease of the mutual information due to the measurement. For the 
critical case a = 0:4, this quantity behaves as ASm ~ l.55d^^'^^ for large d. 



The mutual information represents the total correlation of the bipartite system including 
the quantum and the classical correlations 16|. This quantity is suitable for understanding 
the behavior of QET for the setting 1 because the decrease of the logarithmic negativity 
becomes zero for rf > 1, and it seems that we do not have any quantum correlations between 
A and B. In this case, the information on the entanglement of the ground state is encoded 
as classical correlation between A and B, which is established after the POVM measurement 
of A. As d increases, the mutual information monotonically decreases. For the critical case, 
we have 

ASm ~ l.55d-°-^\ (38) 

The behavior of the mutual information is consistent with the d dependence of the optimized 
energy; the larger the amount of the mutual information, the larger the amount of the 
optimized energy. Thus, we can conclude that the amount of energy teleported via the 
protocol of QET is related to the amount of breaking of the correlation between the group 
A and B. Here it should be stressed that the "classical" correlation supporting energy 
teleportation is induced by ground-state entanglement which is purely quantum. If we have 
no entanglement in the ground state, we do not have any correlation between A and B. 

For the critical case, we obtained the scaling behavior (l37|l and (138|) numerically. However, 
we have not developed any theoretical explanation for the values of these power law indices. 
We expect that an analysis based on conformal field theory may reveal these values. 
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B. Setting 2 

In this setting, the group A consists of A^ — 1 sites and the group B is its complement 
and consists of a single site (Fig. [H]). 



measured sites 

21 + 1 




group A 



group B 



FIG. 6: The setting 2 of our numerical calculation. We change the number of measured sites in 
group A. 

We change the number 2^ + 1 of measured sites in group A. As we explained in Sec. Ill, 
the procedure of the measurement infuses energy into the system. Hence, as £ increases, the 
state of the bipartite system after the measurement is strongly affected by the measurement 
process. In this sense, we can control the strength of the measurement from "weak" to 
"strong" by increasing the number of measured sites i in the group A. We expect that the 
stronger the measurement becomes, the larger the amount of entanglement breaking and 
this enables us to obtain a larger amount of teleported energy of the site B via the protocol 
of QET. 
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< 




FIG. 7: Left panel: decrease of the logarithmic negativity due to the measurement. Right panel: 
optimized energy of site B. 

Figure [7] shows the decrease of the logarithmic negativity due to the measurement, and 
the optimized energy obtained for the site B via QET for the system size N = 100. We have 
also investigated the system sizes 2 < N < 100 and confirmed that the behavior of these 
quantities is same. Hence, we present the result with A^ = 100. As the measurement becomes 
"strong" {£ increases), the amount of entanglement breaking and the resulting optimized 
energy of site B increases. In this setting, as the total bipartite system is pure state, the 
system becomes separable after the measurement when the all sites in A are measured and 
the entanglement breaking due to the measurement is maximal. This behavior is consistent 
with our naive expectation that the ground- state entanglement of the harmonic chain is a 
resource for QET. Figure [8] is the i dependence of the ratio of the optimized energy of site 
B and the entanglement breaking. 
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FIG. 8: Ratio of the optimized energy of site B and the amount of entanglement breaking. 



As we can observe, this ratio is monotonically increasing with respect to the size of the 
measured site I and bounded from above. The maximum value of the ratio is achieved at 
i, = N/2 — 2, which is the maximal number of the measuring sites in group A. Thus, the 
maximal amount of teleported energy is bounded by the amount of entanglement breaking 
and the following inequality holds: 



\{HB)\<f3AE^, (3<1, 



(39) 



where /3 is a constant whose value depends on the system size N and the parameter a. To 
find the relation between the amount of entanglement breaking and the optimized energy, 
we investigated the N dependence of AE^ and {Hb) for a specific value of £ = N/2 — 2; 
this value of i corresponds to the maximal possible number of measuring sites in group A 
for a given A^. The entanglement breaking due to the measurement is expected to become 
maximal (the strongest measurement case). 

As the system size A^ increases, the amount of entanglement breaking and the teleported 
energy both decrease (Fig. [9]). 
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FIG. 9: N dependence of AEn and \{Hb)\ for £ = N/2- 2. 



For the non-critical case a = 01,02,0:3, these quantities approach constant values for large 
A^. For the critical case o = 04, we have the following scaling behavior 



AE 



N 



iN- 



-0.32 



2.1 



(Hb) ~ 0.12N-^-' + 0.0020613. 



By taking the ratio of these quantities, we found that the following inequality holds: 

\{Hb)\</3{N)AE^. 
The behavior of the function (3{N) is shown in Fig. [TUl 
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FIG. 10: The function /3(iV). 

For o = oi, 02, 03, (3 asymptotically approaches constant values. For the critical case o = 04, 
(3 behaves as 

P{N) ~ 0.00026A^°-=^^ (42) 
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The inequality ( BTj) implies that a large amount of teleported energy requires a large amount 
of consumption of the ground- state entanglement between the groups A and B. In other 
words, the possible amount of teleported energy is bounded by the amount of entanglement 
breaking. This relation was previously confirmed for the minimal QET model by one of 
authors ^]. We have confirmed that a similar relation holds for the QET model with a 
harmonic chain. 

It is noted that our bound (HTl) with N = 2 does not quantitatively coincide with the 



3], but a harmonic chain. However, 



bound of [3:] because we treat not a qubit chain as in 

there also exists a qualitative difference between our result with N = 2 and the result in 

Due to noncommutativity between the interaction Hamiltonian and measurement operators, 

our measurement of A disturbs energy of B directly when N = 2. Thus the energy gain 

from B is not purely interpreted as QET output. (However, when A^ > 4, the measurement 

of A does not change energy of B; thus, the energy gain from B relies only on QET.) In |3|, 

the measurement of one qubit is a non-demolition measurement for the other qubit. Thus, 

fl 
even in the case of two qubits, the energy gain of [3| is transported by QET. 

V. SUMMARY 

We numerically investigated the protocol of QET for a harmonic chain model. For a 
bipartite system in the harmonic chain, we applied the QET protocol defined via a POVM 
measurement and LOCC. The resource for QET in the harmonic chain is the ground-state 
entanglement between spatially separated two groups. We can extract energy from the 
system by breaking the entanglement via the POVM measurement and following LOCC. 

We confirmed that the amount of extractable energy via the protocol of QET is bounded 
from above by the amount of the entanglement consumed due to the POVM measurement. 
This implies that a large amount of teleported energy needs a large amount of consumption 
of the ground-state entanglement between the two groups in the bipartite system. We also 
considered the situation that the logarithmic negativity between two groups is zero. In 
this case, the two groups are separable and there is no quantum correlation between them. 
Even in such a case, we showed that we can extract energy via the protocol of QED and 
the amount of energy is correlated with the amount of breaking of the mutual information, 
which quantifies the total correlation between the two groups. It should be emphasized 
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that the above "classical" correlation supporting QET is originally induced from a purely 
quantum correlation, that is, the ground-state entanglement. If we have no entanglement 
in the ground state, such "classical" bipartite correlation between A and B does not take 
place at all. As pointed out above, QET only needs "classical" correlation of A and B. This 
result suggests that QET processes are tolerant of decoherence which destroys the quantum 
correlation of A and B, unlike the usual quantum teleportation. 
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